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I present an extension to arbitrary A^ of a previously proposed field theoretic 
model, in which unitary amplitudes for 1 ^ 8 processes were obtained. The Born 
amplitude in this extension has the behavior A^^^j^ = g N\ expected in a 



^ ■ bosonic field theory. Unitarity is violated when |yli_>7v| > 1, or when A^ > Acrit — 

e/g. Numerical solutions of the coupled Schrodinger equations shows that for weak 

^ ■ coupling and a large range of A^ > Acrit, the exact unitary amplitude is reasonably 

H ' 

5^ ! fit by a factorized expression 

|Ai->iv|max ^ 0.73 • - ■ exp (-0.025/(72) . 

The very small size of the coefficient l/g"^ , indicative of a very weak exponen- 
tial suppression, is not in accord with standard discussions based on saddle point 
analysis, which give a coefficient of ~ 1. The weak dependence on A^ could have 
experimental implications in theories where the exponential suppression is weak 
(as in this model). Non-perturbative contributions to few-point correlation func- 
tions in this theory would arise at order K ^ ((0.05/(7^) -f 2 InN^ / ln{l/ g"^) in 
an expansion in powers of g"^. 



Tree-level calculations of amplitudes for scattering processes in which many 
bosons are produced fail to obey unitarity when the number produced is too large. 
Specifically, in either massive scalar [1,2] or vector [3] field theories, the 1 ^ A^ 
tree amplitude for fixed angle scattering has the unacceptable behavior 

for large A^, where g'^ is the quartic coupling. Cross sections then violate unitar- 
ity when A^ ^ l/fi*^ for E^N large, E/N finite. As shown explicitly in refs.[2,3], 
the problem originates in the coherence of the approximately A^! graphs which 
contribute. The cure surely lies in the summation of all loops, since their contri- 
bution is equal to that of the tree graphs precisely when the latter become large. 
The question remains: does the unitary damping suppress the cross section at an 
exponentially small value? or do coherence effects allow an experimentally inter- 
esting cross section for producing many massive bosons in a high energy collision? 
Independent of the phenomenological application lies the question as to how to 
calculate high energy, multiparticle processes in field theory. 

Zakharov [4] has given a thoughtful argument for exponential suppression, 
which I will simply paraphrase. Consider the dispersion relation for the fourier 
transform of some two point function in field theory or quantum mechanics, with 
enough subtractions to make it convergent. At zero energy, we have 

m) = ^I^ImU{E) , (1) 

m 

where 

Imn(E) -^ lAi^iVp^(^-miV) • (2) 

iV 

Suppose we now assume that n(0) has an asymptotic expansion in coupling 

n(o) = J2^k g'' , (3) 



where for large K 

K 
n(0) - J] afe ^2fe 



K 



< {g'/hr K\ . 
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Then, by the usual arguments, the error in truncating the series is minimized at 
K = Kcrit = b/g'^, and the error in omitting the remainder is bounded by 

An(0) < e-^/3' . (4) 

At this point the discussion may take the following form: 

Suppose we either know or think we know the value of b. (For example, in the 
scalar (p^ field theory, it is often assumed that b is given correctly as IGvr^ 
by the saddle point analysis of Lipatov [5].) Then the high order behavior 
(for k > b/ g^) of Im 11 is bounded by exp (—b/g"^). Assuming that the high 
multiplicity n > b/g piece of Im 11 is an important component of the high 
order piece, one concludes that high multiplicity contributions to Im 11 (or 
the cross section, in field theories) are bounded by [6] 

E 1^1-^1' < ^"'/'^ • (5) 

N ^ b/g'2 

One may also turn the argument around: suppose that we have some evidence 
that multiparticle contributions to Im 11 are as large as exp (— c/(?^). Then, from 
the dispersion relation, we may conclude that the perturbation series for 11(0) 
becomes untrustworthy for k ^ c/g"^. If c ^ 1, the situation becomes interesting, 
on two counts : (1) although formally exponentially suppressed, there may be 
experimentally interesting multiparticle cross sections if g"^ is not too small (2) as a 
consequence of the dispersion relation, the two point function may have interesting 
non-perturbative contributions at relatively low order k ~ c/g . 

In this paper, I will propose an extension to arbitrary A^ of a field theoretic 
model [7] previously applied to a calculation of a unitary 1^8 transition ampli- 
tude. In this extension, the (time-dependent) amplitudes Ai^jv iiiay be directly 



(numerically) calculated . The result of the present calculation is that in the weak 
coupling limit, the amplitude Ai_,jv is very weakly exponentially suppressed: 



A,^N - e-O-025/.^ , (6) 

with the implications mentioned in the preceding paragraph and in the abstract. 

A field theoretic model, loosely based on a field theory, has been presented 
[7] as a laboratory for the study of multiparticle amplitudes. The model is defined 
through its Hamiltonian 

oo oo 

^ = X^ °I°fc + 55' ^ $Z 4+fe ^^ (^k'P + h.c. . (7) 

fc=l j,k=l 

The modes labeled by i,j, k will be called "momenta", and the action of the her- 
mitean projection operator P on a state vector {ip) is as follows: 

P 1^) = 

if there is more than one quantum in the state with any given momentum and 



v\ 



otherwise. V has been introduced in order to approximate the infinite volume 
effect of field theory (in box normalization): namely, one generally omits consid- 
eration of amplitudes for transitions to states with more than one particle in a 
given (discrete) momentum state. Thus, we exclude "laser" effects. In this sense, 
there is an exclusion principle without imposing anticommutation relations and 
antisymmetrization. Other than that, H resembles a 0"^ field theory in a cavity, 
with no temporal or spatial derivatives in the lagrangian. It is also a kind of matrix 
model. 



It is an important consequence of (7) that the momentum operator 



fc=l 



k a\ak 



is a constant of the motion. Thus, the Hilbert space factorizes into subspaces 
with definite P. Because of the positivity of all of the momenta, these will he finite 
dimensional subspaces. Within a subspace of given P, there will be sectors charac- 
terized by different numbers of particles n. A subspace with definite P may contain 
a maximal state, namely an A^-particle state with momenta k = 1,2, ... N, such 
that P = N{N + l)/2. I will restrict my study to these subspaces, which can be 
labeled by A^. For A^ = 8, the number of states with n = 1,2, .. .8 particles is 
(1, 17, 91, 206, 221, 110, 21, 1) respectively. A numerical analysis of this case [7] 
proceeded as follows: with the system initially in the one-particle state, the 668 
coupled (complex) time- dependent Schrodinger equations were numerically inte- 
grated in order to obtain the maximum 8-particle amplitude |A8(t)|max- The result 
indicated an approximate behavior for the 1^8 transition 

|^T-8|max — c 

for g'^ > 0.07. This is the value of g"^ for which the Born amplitude violates unitar- 
ity: |Af°™(t)|n,ax>l. 

In general, the total number of states for a given P can be obtained as the 
exponent of x in the expansion of the generating function 

00 CX) 

l[{l + x^) = Y,Afpx'' . (9) 

i=l P=0 

This reveals that even an extension to A^ = 9 means quadrupling the subspace from 
668 to 2048 states. Thus, exact numerical analysis for larger A^ rapidly becomes 
impractical. 



How then to extend the results beyond N = 81 Clearly, an approximation to 
the dynamics of the model at large A^ is required; however, it is imperative that 
the extension to arbitrary A^ (a) retain the unitarity property and (6) display the 
gN-l jYj behavior in the tree, or Born approximation. 

A detailed study of the coupling structure of the Hamiltonian (7) for A^ = 8 
yields a strong clue as to how this extension may be accomplished. If one asks, 
how many states (on the average) with n — 1 particles does a particular state with 
n particles couple to, then one arrives at the following list: 



8- 


-^7 


16 states 


7- 


-^6 


15± 1 states 


6- 


-^5 


12± 1 states 


5- 


-^4 


9± 1 states 


4- 


-^3 


6 states 


3- 


-^2 


3 states 


2- 


■^ 1 


1 state 



If we ignore the 8^7 result, and are liberal with the 7^6, then it is not too 
far amiss to conclude that, approximately, each state with n particles couples to 
an average of n{n — l)/2 states with n — 1 particles. [8] This appeals to intuition: 
adding each (unequal) pair of 'momenta' in the given state with n particles will, 
in large probability, give one of the possible states with n — 1 particles. Thus, I 
propose as an extension to (7) a quantum mechanical system in which there are 
simply A^ states labelled by n = 1, . . . A^, differing in energy by a constant amount 
(taken to be 1 in our units), and in which each state \n) is coupled with strength 
g n{n — l)/2 to the state |n — 1) . In the interaction representation the dynamics 
of this model is embodied in the system of time dependent Schrodinger equations 

^An{t) =g [^^^^ A,,_i(t) e^* + ^^^^ A„+i(t) e^^ 

(10) 
n = l...N 

where Ao(t) = A7v+i(t) = 0. For a given A^, these equations may be solved by first 



diagonalizing the N xN Hamiltonian, or by direct numerical integration. But first, 
I state without details of proof a result which establishes this system as a viable 
laboratory in which to study the large- A^ problem: 

If An{0) = 6ni, then in Born approximation 



iBorn/j.N „n—l ' ^ ^ 



AT''{t)=g--' (^-- n! . (ll^ 



This can be proven by direct subsitution into Eq.(lO) with no rescattering {i.e., 
dropping the second term on the right). It was actually arrived at constructively 
by using Laplace transforms. 

As a result of (11), 

l^r^lmax = g"""' N\ . (12) 

I note that in the exact version of the model, with N = 8, [7] it was found that 

lAf-Uax ^0.27/8! . (13) 

Note also that in Born approximation the anharmonic oscillator with coupling 
jg'^x'^ gives [9,10] 

iBorni „N~1, 



\^N Imax ~ g ViV! . (i4j 

In both the present case and in the case of the anharmonic oscillator the state |A^) 
with unperturbed energy A^ is normalized to unity. 

Since there is a normalization condition ^^ |An(t)| = 1, it is seen from Eq.(12) 
that unitarity will certainly be violated in Born approximation when 

N> Merit = - . (15) 

g 

For a given g , this defines the value of A^ for which rescattering terms are essential. 



It is then a simple matter to go back to Eq. (10), and find (numerically) the 
maximum value |Ajv|max attained by Aj^(t) as a function of N, g^. [11] 

Results : The results of this work are encompassed in Figs. (1) and (2). In Fig. 
(1), I display on a log — log plot the A^ behavior of |Ajv|max, for a large range of 
values of g"^. The range of A^ for each g"^ is chosen to lie above A'^c^it. Two striking 
observations may be made from the graph: (1) The A^-behavior is universal over 
the whole array of g — the amplitude factorizes; (2) The curves are excellently 
fit with a simple inverse proportionality I^Tvlmax ~ 1/A^. The wavy oscillations 
in the curves for the larger values of g"^ are real: they can be seen in detail in Fig. 
(3) as an odd-even effect. 

In Fig. (2), I plot I^Tvlmax vs. Xj g^ on a log plot for a range of values of A^. 
Again, one notes the factorizability. For the smallest values of (7^, the curves are 
well fit by the exponential form |Ajv|max ~ exp (— 0.025/(?^). In sum, therefore, 
the results of the numerical study of the equations (10) is that 

l^i^ivlmax ^ 0.73 • ^ ■ exp (-0.025/(?2) . (16) 

I close with some discussion of these results. 

Remarks and Conclusions: (1) First, the factorizability. The results of Lipatov 
[5] and the graphical analysis of Parisi [12] support the notion that the contribution 
of a high order K of perturbation theory to an A^-legged Green's function should 
be (roughly) independent of A^, for A^ not too large. What I find here is that the 
functional dependence on coupling constant of the exact non-perturbative A^-point 
function is independent of A^, in the region where A^ is large (A^ > e/^f.) 

(2) The behavior of Ai^tv in the model being examined is totally different from 
its behavior in the case of the anharmonic oscillator. In that case, for Ng'^ not too 
large, I Aat I max ~ exp(-A^). [9,10] 



(3) Another surprising result is the generation of a small dimensionless number, 
namely 0.025, the coefficient of 1/ g^ in the exponential suppression. Typical saddle 
point analyses (even in the theory [13]) give coefficients of 0(1) {modulo factors 
of (47r) for d dimensions). The implication is that for A^ > Acrit, the theory 
becomes strongly coupled for very small values of g'^. The next point is related to 
this one. 

(4) The contribution of the 1 — * A^ excitation to a few-point function (such as 
the 1 — * 1 transition amplitude) can be estimated as 

(1, t = T\l, t = 0)^ ~ \An\L. . (17) 

If we wish to see at what order K the non-perturbative contribution of the right- 
handed side competes with a perturbative development of the left-hand side, we 
can use the parameterization (16) and set 



27^cHt ^ exp (-0.05/^2)/;y2 ^ ^^^ 
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with A^ > e/g. This gives the formula in the abstract. Even for g'^ = 0.005 and 
A^ = 100, this gives Kent = 3.6, a very low value. 

(3) Finally, a comment on the limitations of the numerical study. As A^ be- 
comes large and/or g"^ becomes small, the roundoff errors become more impor- 
tant. At some point, |^iv|max is smaller than the deviation from unitarity 5 = 
|^j^|Aj^(t)p — ll caused by finite numerical accuracy. This occurred for g'^ < 
0.005 (and A^ > Acrit)- The results for such small values of coupling may therefore 
be untrustworthy, and at present I am examining alternate methods of exploring 
the full range of g in the non-perturbative region. 

To conclude, I have found in a unitary model a large 1 -^ A^ amplitude which 
is only very weakly exponentially suppressed. The dependence on A^ is also very 
weak. This may hold out some possibility for observing large multiplicity central 
events in high energy collisions. 
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Figure Captions 

Fig. 1 The maximum value attained by |A7v(t)| as a function of A^, for various 
values of g^. The Born approximation fails for A^ > A^crit = ^/ 9- For small 
(^2, large N, |Aiv|max ~ l/iV- 

Fig. 2 |^i-+iv|max vs. 1/ g'^ for three values of A^. The behavior in 1/ g'^ is inde- 
pendent of A^ : Ai^N ~ exp (-0.025/5-2). 

Fig. 3 Detail of Fig. (l) ior g"^ = 0.10, showing origin of waviness as an odd-even 
effect. 
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